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Abstract 



Non-leading contribution to the pion electromagnetic form factor which comes from 
the pion twist-3 wave function is analyzed in the modified hard scattering approach 
(MHSA) proposed by Li and Sterman. This contribution is enhanced significantly 
due to bound state effect (the twist-3 wave function is independent of the fractional 
momentum carried by the parton and has a large factor ~ m^/mo with being 
the pion meson mass and uiq being the mean u- and d-quark masses). Consequently 
although it is suppressed by the factor l/Q 2 , the twist-3 contribution is comparable 
with and even larger than the leading twist (twist-2) contribution at intermediate 
energy region of Q 2 being 2 ~ 40GeV 2 . 



(N 
> 
m 
O 

(N 



On 




X 



1 Mailing address. E-mail: caofg@itp.ac.cn. 



1 



1 Introduction 

There has been a lot of discussions about applying perturbative QCD (pQCD) to exclusive 
processes at large momentum transfer [1-15]. Although there is general agreement that 
pQCD is able to make successful predictions for the exclusive processes at asymptotic limit 
(Q 2 — ► oo), the applicability of pQCD to these processes at experimentally available Q 2 
region has been being debated and attracted much of attention. The difficulties in practical 
calculation mainly come from the end-point singularity, i.e. in the end-point region (x — *■ 0, 1 
with x being the fractional momentum carried by the parton) the virtuality of intermediate 
states is small and the running couple constant a s becomes large, thereby perturbation 
expansion might be illegal. However, perturbative calculation can be rescued with the help 
of some techniques to cure the end-point singularity [8-15], for example, the incorporation of 
the transverse structure of the pion wave function || |9|, pL0[| , the introduction of an effective 
gluon mass |ll[ and a frozen running coupling constant [JTT], |I2| . Recently, Li and Sterman 



fT3| , |l4f proposed a modified hard scattering approach (MHSA) for the hadronic form factor 
by taking into account the customarily neglected partonic transverse momentum as well as 
Sudakov corrections. They point out that pQCD calculation for the pion form factor begins 
to be self-consistent at about Q ~ 20Aqcd, which is similar to the conclusion given in Ref. 
|| . More recently, Ji, Pang and Szczepaniak [|l5j arrived at a similar conclusion as Refs. 



13] , [14| by analyzing the factorization perturbation formalism for the pion form factor 



in the framework of light-cone time-order perturbative theory. These studies shed light on 
applying pQCD to exclusive processes at intermediate energy region. 

However, there is still a crucial problem which has not been solved, that is although 
improved pQCD calculation for the exclusive processes is self-consistent at currently exper- 
imentally accessible energy region, the numerical predictions are generally far smaller than 
the experimental data. For example, pQCD prediction for the pion form factor is 



i^(Q 2 ^oo) = 167ra s (Q 2 )CV J[dx][dy]<l>(x)^-^<f>{y) 
167ra s (Q 2 )f 2 



Q 2 



where [dx] = dxidx2S(l — x\ — x 2 ), [dy] = dyidy26(l — yi — y 2 ), f n = 93 MeV is the pion decay 
constant, and <f)(x) is the distribution amplitude of the pion meson. The asymptotic form 
for the distribution amplitude has been employed in obtaining the send expression in Eq. 
(H), since any distribution amplitudes for the pion meson should approach the asymptotic 
form as Q 2 — ► oo, 
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X 



v / 3/ 7r x 1 x 2 . (2) 



Eq. (H) gives only 1/3 of the experimental data at intermediate energy region. Although 
the Chernyak and Zhitnitsky (CZ) model for the distribution amplitude 

^ cz \ x ) = 5V3f 7T x 1 x 2 (x 1 - x 2 ) 2 (3) 
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may enhance the prediction for the pion form factor to the correct direction, the perturbative 
calculation with CZ distribution amplitude has been criticized seriously || [7| because the 
nonperturbative end-point region is much emphasized in the CZ model. Recently studies on 



the pion-photon transition form factor [16| also show that the pion distribution amplitude 
at currently experimentally available energy region is much like the asymptotic form but not 
the CZ form. Hence, how to match the perturbative calculation with the experimental data 
is an interesting issue. There are two possible explanations: one is that non-perturbative 
contributions will dominate in this region; the other is that non-leading order contributions 
in perturbative expansions may be also important in this region. To make choice between the 
two possible explanations one needs to analyze all of the important non-leading contributions 
carefully. These contributions come from higher-twist effects, higher order in a s and higher 
Fock states etc.. Field, Gupta, Otto and Chang [IT] pointed out that for the pion form 
factor the contribution from the next-leading order in a s is about 20% ~ 30%. Employing 
the modified hard scattering approach []13~ , |14}| , Refs. |18[ and [[0| considered the transverse 
momentum effect in the wave function and found that the transverse momentum dependence 
in the wave function plays the role to suppress perturbative prediction. More recently, Tung 
and Li [^U] reexamine the perturbative calculation for the pion form factor in the MHSA by 
respecting the evolution of the pion wave function in b (the transverse extent of the pion) and 
employing the two-loop running coupling constant in the Sudakov form factor. It is found 
|20| that the evolution of the pion wave function in b improves the match of perturbative 
prediction with the experimental data. However, in order to answer the question whether 
the perturbative calculation is able to make reliable prediction for the exclusive processes at 
currently experimentally available energy region, the other non-leading contributions such as 
that from higher twist effects and higher Fock states |2|, 21 should also be analyzed carefully. 

It has been expected that the power corrections to the pion form factor (~ 1/Q 4 ) which 
come form the higher twist terms of the pion wave function may be important in the inter- 
mediate energy region [|22], [23], |2|, |5|, since there is a large factor ~ / m (m w being 
the pion meson mass and m being the mean u- and d-quark masses) in the twist-3 wave 
function. However, the calculations for these higher twist contributions are more difficult 
than that for the leading twist (twist-2) because of the end-point singularity becoming more 
serious. The leading twist wave functions in the initial and final states being proportional 
to X1X2 (xi and X2 being the fractional momenta carried by the quark and anti-quark) 
and 2/12/2 (see Eq. (fj)) may cancel the end-point divergent factor I/222/2 coming from the 
hard-scattering amplitude. However, the asymptotic behavior of twist-3 wave function is x- 
(y-) independent (see Eq. (^)), which has no help at all to cure the end-point singularity. 
In this case, Sudakov form factor is expected to be able to assure the reasonableness of the 
perturbative calculation. Unfortunately, the estimations for the twist-3 contribution in the 
medium energy region do not agree with each other [22, |2^, 24, 25. |2(| (see Fig. 1). Ref. 
predicates that 
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with 



J(Q 2 



lnln(Q 2 /A| CD ) + 



(5) 



The first and second terms in Eq. (f|) correspond to the leading twist (twist-2) and next-to- 



leading twist (twist-3) contributions respectively. In Ref. [^3|] the double logarithmic (DL) 
corrections are calculated in the one loop approximation and it is supposed that the sum of 
all DL corrections transforms to the exponential function form (Sudakov form factor). Hence 
it is argued that the divergent factor \jx 2 y 2 at £2(3/2) - is modified by the following way 



X2V2 



— exp {-[a s (Q 2 )/27r]C F L(x 2 , y 2 , k ± , Lj} 



X2V2 



(6) 



with 



L(x 2 ,y 2 ,k ± ,l ± ) = ln(Q 2 /kl)ln(l/x 2 ) -ln 2 (l/x 2 ) + ln(Q 2 /ll)ln(l/y 2 ) -ln 2 (l/y 2 
-ln(l/x 2 )ln(l/y 2 ). 



(7) 



It is argued [p3j that the integral with function L(x 2} V2,k± } l±) can not be calculated unam- 
biguously. This uncertainty is incorporated to the factor a being 1 < a < 2 in the function 
J{Q 2 ) (Eq. (H)). According to Eq. (^), the twist-3 contribution is larger than the asymptotic 
term (twist-2 contribution) in the region of Q 2 < 30 GeV 2 . Ref. includes the Sudakov 
corrections in a similar way as Ref. [23] but improved the estimation on the function J(Q 2 ) 
and the running mass m(Q 2 ), and gives 



167ra s (Q 2 )f 2 

Q 2 
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(8) 



It can be found from Eq. (§) that the twist-3 contribution is larger that the twist-2 contri- 
bution at about Q 2 < 15 GeV 2 . Ref. [[25j analyzes the Sudakov effects by introducing an 
cut-off on the integral region instead of introducing the transverse momenta k± and l±, and 
gives another prediction 



167m s (Q 2 )/ 2 

Q 2 




Q 2 m\ 6 
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(9) 



Eq. (H) tells us that the twist-3 contribution is about 2 ~ 0.6 of the leading twist contribution 
at the energy region 2 GeV 2 < Q 2 < 10 GeV 2 . All of the above calculations (Eqs. (§, (|) 
and (§)) give correct power suppression (~ 1/Q 2 ) behavior for the twist-3 contribution in 
the large Q 2 region, but their predictions for the dependence on InQ 2 are very different. 
The main reason for these differences is that Sudakov corrections are evaluated in different 
approximations in Refs. f53]| , |f24|j and | 2"5| . In the modified hard scattering approach for the 
exclusive processes proposed by Li and Sterman |L3], [HJ] , the customarily neglected partonic 
transverse momentum are combined with Sudakov corrections, and the Sudakov form factor 
is expressed in a more convenient space (6-space), which provides an more reliable and 
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systematical way to evaluate the Sudakov effect. Li and Sterman's formalism is originally 
obtained for studying the contribution from the leading twist wave function. We point 
out that for the pion electromagnetic form factor the MHSA can be extended to evaluate 
the contribution coming from the twist-3 terms of the pion wave function. One manifest 
advantage of MHSA is that there is no other phenomenological parameter but the input 
wave function need to be adjusted. The purpose of this work is to analyze the twist-3 wave 
function contribution to the pion form factor in the framework of MHSA. 



2 Formalism 



We first review the derivation of the modified hard-scattering formalism for the leading 
twist (twist-2) contribution to the pion form factor Taking into account the transverse 
momenta and lj_ that flow from the wave functions through the hard scattering leads 
to a factorization form with two wave functions ip(x, kj_) and ijj(y, corresponding to the 
external pions combined with a hard-scattering function T H (x, y, Q, k^, which depends 
in general on transverse as well as longitudinal momenta, 

Ft 2 \Q 2 ) = J[dx}[dy] J rf 2 k ± rf 2 l^(^k ± ,P 1 , / x)Ti* =2) (x,y,g,k ± ,l ± ,/i)^(i/,U,P 2 ,/i), (10) 

where P± and P 2 are the momenta of the incoming and outgoing pion respectively, Q 2 = 2 Pi ■ 
P 2 and fj, is the renormalization and factorization scale. To the lowest order in perturbation 

(t=2) 

theory, the hard-scattering amplitude Tjj is to be calculated from one-gluon-exchange 
diagrams. Neglecting the transverse momentum dependence in the numerator of T^ -2 ^ one 
can obtain, 



Th 2 \x,y,Q,k ± ,\ ± ,fi) 



x 2 Q 2 



[x 2 y 2 Q 2 + (kj 
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where Cp = 4/3 is the color factor and a s (fi) is the QCD running coupling constant. The 
first and the second terms in the denominator come from fermion and gluon propagators 
respectively. 

Eq. ( |10|) can be expressed in the b- and h-configurations via Fourier transformation 



Ft 2) (Q 2 



[ dx \ [dy] , n no, ^ 9 <p(x, b, Pi,fi)T^ =2 \x, y, Q, b, h, fi)(p(y, h, P 2 , fi), (12) 



(2tt) 2 (2tt) 2 ' 

where wave functions <p(x, b, Pi, fi) and (p(y, h, P 2 , //) take into account an infinite summation 
of higher-order effects associated with the elastic scattering of the valence partons, which 



give Sudakov suppressions to the large-fo(/t) and smal\-x(y) regions |T3], |27], ^8 

-s(Z,b,Q)-s(l-Z,b,Q)-2 [ 

Ji 



y(f,b,P,/x) = exp 



n/b [i 



x 



(13) 
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-a s /iT is the quark anomalous dimension. s(£,b,Q) is Sudakov exponent factor 
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ln 2 (2g) - ln 2 (-26) 



(14) 



where 
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rif is the number of quark flavors and 7 is the Euler constant. In the derivation of Eq. 
(fUj), the one-loop running coupling constant has been employed. It is pointed out |29[ that 
additional two terms will appear in the s(£, 6, Q) expression if the two-loop running coupling 
constant is used. The two terms reduce the prediction for the pion form factor by only a few 
percent p9| in the intermediate energy region. So for simplicity, we neglect these terms. 

Applying the renormalization group equation to T# and substituting the explicit ex- 
pression for Tff , one can obtain the following expression for the pion form factor 

Ft 2) (Q 2 ) = 



[dx][dy] j bdb J hdh \QT{CF(^s{t)x2Q 2 K {^/^Qh)(p{x,l/b)(p{y,l/h) 
[6(b - h)K (^Qb)I (^Qh) + 6(h - b)I (V^Qb)K (^Qh)] 

(16) 



x exp(-S(x,y,Q,b, h,t)), 



where 



S(x,y,Q,b, h,t) 



E 

\i=l 



S(Xi 



b,Q)+i2s(Vi,Q,h)) - 1 In - i- In -t 
i= i J Pi -b Pi —h. 



(17) 



K and Iq are the modified Bessel functions of order zero, t is the largest mass scale appearing 
in T 



(t=2) 
H 



t = max {^JxyQ, 1/6, 1/h) . 



(18) 

If b is small, radiative corrections will be small regardless of the values of x because of the 
small a s . When b is large and xyQ 2 is small, radiative corrections are still large in T^~ 2 \ 
but ip will suppress these regions. In Eq. fll6|) , (f>(x, 1/6) and (f)(y, 1/h) are two input "wave 
functions" which respect the non-perturbative physics. In the large-Q 2 region, they can be 
taken as the asymptotic form of the twist- 2 distribution amplitude (Eq. (fj)) |13|, |27|, |28fl . 

In the above discussion, only the leading twist wave function is considered. Now, we 
address the contributions coming form the twist-3 wave functions. The operators which 
contribute to the twist-3 parts of the pion wave function include 75 and 7507^, and the two 
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matrixes might mix under the consideration of the evolution equation for two-quark state in 



the pseudoscalar channel. It is pointed out in Refs. £24| that the twist-3 wave function 
of pion can be expressed as 



■0 



(*=3), 



X, k ± ) ^ 75<i 



l-i 



(19) 



where is the partonic transverse momentum. 3 is the distribution amplitude of twist-3 

where = 139 MeV is the pion meson mass and m(Q 2 ) is the mean value of the u- and 
d-quarks masses at the scale Q 2 , 



m{Q 2 



' a s (Q 2 ) \ 



2\ \ 4//3o 



(21) 



with /3b = 11 - §n/, and m (l GeV 2 ) = 7 ± 2 MeV. 

The hard scattering amplitude for the twist-3 wave function differs from that for the 
twist-2 wave function, which turns out to be 



Th 3) (x,y,Q,k ± ,l±,n) 



64iiCf d s {jj)x2 



x 2 Q 2 + k ± 2 [x 2 y 2 Q 2 + (kx - lx) 2 ] 



(22) 



Following the derivation for the leading twist wave function we can obtain the twist-3 con- 
tribution to the pion form factor in the modified hard-scattering approach, 



F^\Q 2 ) = J[dx][dy] I bdb J hdh647rC F a s {t)x 2 K {^ 2 -Qh)(f> 3 (x)(f> 3 {y)exp(-S(x,y,Q,b,t)) 
x [e{b-h)K {^ 2 Qb)h{^ 2 Qh) + e{h-b)I,{^ 2 Qb)K a {^y- 2 Qh)}. (23) 



It can be found that the hard scattering amplitudes T^ -2 "* (Eq. [TTJ) and T^~ 6> (Eq. |22[) are 
divergent in the end-point region x, y — *■ 0, k±,l± — > 0. However, the twist-2 contribution 
to the pion form factor can be calculated readily because the twist-2 wave functions being 
proportional to X\X 2 and y\y 2 respectively cancel the divergent factor l/x 2 y 2 in the T^~ 2 \ 
Furthermore, the Sudakov corrections also suppress the contribution from the end point 
region. For the twist-3 contribution, the wave function is a constant in the whole region of 
x (see Eq. fl20D ), which has no help at all to cure the divergent factor l/x 2 y 2 in the T^~ 3 \ 
In this case, the Sudakov form factor may guarantee that the calculation is reliable since the 
factor e~ s rapidly decreases to zero more rapidly than any power of x(y) at the end-point 
region (see Eqs. ( fH| ) and (|TTD ) . 



<t=3) 



3 Numerical result and discussion 

We present the numerical evaluations for the twist-2 and twist-3 contributions to the pion 
form factor in Fig. 1. The thinner solid curve is MHSA prediction for the twist-2 contribution 
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(Eq. fllBD). The thicker solid curve is twist-3 contribution to the pion form factor obtained 
in this work (Eq. (|23|) ), while the dash-dotted is the result of Ref. p3| (the second term 
in Eq. (^) with a = 1.5). The dotted and dashed curves are the results given in Refs. p4| 
(the second term in Eq. (||)) and |25| (the second term in Eq. @) respectively. All of 
the calculations given in this work, Refs. [p3|| , p4[ and |25| show that compared with the 



leading twist contribution, the twist-3 contributions are suppressed by the factor 1/Q 2 at 
asymptotic limit (Q 2 — > oo). But the predictions are different in the medium and lower 



energy regions. Our result is much larger than the result of Ref. [[25[] in the energy region 



of 2 GeV 2 < Q 2 < 40 GeV 2 , and a little larger than the result of Ref. || as Q 2 > 5 GeV. 



12 ^ r 



Ref. p3f and this work give very similar results in the large-Q 2 region, but our prediction is 
a little smaller than the result of Ref. [ 23 1 at about Q 2 < 15 GeV 2 . The Sudakov corrections 



are respected systematically in MHSA, while they are evaluated in various approximations 
in Refs. p3[ |24], p5| , so the prediction in this work is more reliable. In Fig. 2, we include 
both twist-2 and twist-3 contributions (obtained in this work) to the pion form factor, and 
compare with the experimental data. The dotted and dashed curves are twist-2 and twist-3 
contributions respectively, and the solid curve is the sum. Compared to the leading (twist-2) 
contribution, the twist-3 contribution is negligible at asymptotic limit since it is suppressed 
by the factor 1/Q 2 . However, the twist-3 contribution is comparable with and even larger 
than the leading twist contribution at intermediate region of Q 2 (2 GeV 2 < Q 2 < 40 GeV 2 ). 
Also it can be found that the perturbative calculations including both twist-2 and twist-3 
contributions are larger than the experiment data at about Q 2 < 5 GeV 2 . One can expects 
that the other nonleading contributions such as those coming form higher Fock states may 
be also important at lower energy regions. 

In summary, we analyzed the twist-3 contribution to the pion electromagnetic form fac- 
tor in the modified hard scattering approach in which Sudakov corrections are respected 
systematically, and compared with various approximate calculations. It is found that the 
twist-3 contribution is enhanced significantly since the twist-3 wave function is independent 
of the fractional momentum carried by the parton and has a large factor ~ m 2 /m , while 
the twist-2 wave function is proportional to X\X 2 (2/12/2) which cancels the end-point diver- 
gent factor I/X2I/2 in the hard-scattering amplitude. Thus, although it is suppressed by the 
factor 1/Q 2 as compared with the leading (twist-2) contribution, the twist-3 contribution is 
comparable with and even large than the leading twist contribution at intermediate region 
of Q 2 being 2 ~ 40 GeV 2 . The perturbative predictions including both twist-2 and twist-3 
contributions are larger than the experiment data at lower energy regions, which indicates 
the importance to study the other nonleading corrections at these energy regions. 
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Figure caption 

Fig. 1 Twist-2 and twist-3 contributions to the pion form factor. Each curve is explained 
in the text. 

Fig. 2 Perturbative prediction for the pion form factor including both twist-2 (dotted curve) 
and twist-3 (dashed curve) contributions. The solid curve is the sum of twist-2 and 
twist-3 contributions. The data are taken from Ref. []5C[] . 
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